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We present a detailed analysis of the nature of electronic eigenfunctions in 
one-dimensional quasi-periodic chains based on a clustering idea recently in- 
troduced by us [Sil et al., Phys. Rev. B 48, 4192 (1993) ], within the frame- 
work of the real-space renormalization group approach. It is shown that even 
in the absence of translational invariance, extended states arise in a class of 
such lattices if they possess a certain local correlation among the constituent 
atoms. We have applied these ideas to the quasi-periodic period-doubling 
chain, whose spectrum is found to exhibit a rich variety of behaviour, in- 
cluding a cross-over from critical to an extended regime, as a function of the 
hamiltonian parameters. Contrary to prevailing ideas, the period-doubling 
lattice is shown to support an infinity of extended states, even though the 
polynomial invariant associated with the trace map is non- vanishing. Results 
are presented for different parameter regimes, yielding both periodic as well 
as non-periodic eigenfunctions. We have also extended the present theory 
to a multi-band model arising from a quasi-periodically arranged array of 
5-function potentials on the atomic sites. Finally, we present a multifractal 
analysis of these wavefunctions following the method of Godreche and Luck 
[ C. Godreche and J. M. Luck, J. Phys. A Math. Gen. 23, 3769 (1990)] to 
confirm their extended character. 
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I. Introduction 

In the last decade there has been intense activity in theoretical studies of 
aperiodic lattices in one dimension. Apart from their intrinsic interest, these 
studies have been mainly motivated by the remarkable discovery of the so- 
called quasicrystalline ordering in solids, beginning with the work of Shecht- 
man et al. 1 in 1984. On the experimental side, one dimensional aperiodic 
lattices have also become interesting objects of study, following the success- 
ful fabrication of layered semiconducting superstructures grown epitaxially 
in accordance with the rules of the well-known Thue-Morse and Fibonacci 
sequences 2 . Parallel with this development, theoretical studies have been 
made on a number of such one-dimensional systems for phonon, electron and 
magnon spectra 3 . While most of these studies have been numerical, there 
have also been quite a few analytical attempts, starting with the pioneering 
work of Kohmoto et al. 4 on the unusual electronic properties of the Fibonacci 
lattice. 

The behaviour of the electronic eigenstates in one-dimensional lattices 
are well-known in the two extreme limits of random and perfect periodic 
ordering; in the former case all the eigenstates are localised, while in the latter 
they are extended. It is precisely because aperiodicity is a level of ordering 
intermediate between these two extreme limits that the electronic eigenstates 
in aperiodic linear systems often are found to be critical 4 . Interestingly, there 
is also evidence for extended states in some of these systems 5-9 . While the 
reason for the existence of critical states may be assigned to the lack of 
periodicity in these systems, extended states arise due to other reasons, as 
discussed below. 

In a recent work 5 we analysed this problem in one-dimensional systems 
and found that the basic reason may be traced to the existence of a particular 
type of correlation among the atoms in the lattice. For the sake of clarity 
we summarise here the main ideas developed in that work. In an elementary 
problem of Anderson localisation of the random distribution of A atoms on a 
host lattice of B atoms, it is well-known that there is an extended eigenstate 
at an energy €a in this system if the A atoms always occur in pairs, that is 
as dimers. This is an instance of a definite kind of correlation leading to an 
extended state in a disordered system 6 . In this case the correlation consists 
in the fact that an A atom always appears as a member of a dimer. We 
may generalise this idea as follows. Consider the distribution of a certain 
well-defined cluster of atoms on a host lattice in one dimension. Suppose 



that this cluster consists of a finite number of identical building blocks , each 
of which is a smaller unit consisting, in general, of several atoms. In the 
above example, for instance, a dimer may be regarded as a cluster of two 
consecutive A atoms. If the building blocks are to be found only within the 
clusters and nowhere else on the lattice, then it can be shown that for certain 
energy values the contribution to the total transfer matrix for the chain from 
these clusters is unity. Therefore, at these special energies only the host 
lattice is effective in determining the electronic properties. If the host lattice 
is periodic, then there is a possibility of extended states at these energies for 
the whole sytem. 

The special energies for which the states are extended are a property of the 
cluster only and do not depend on the lattice as a whole. Thus in our earlier 
example of random AA dimers on a B type host lattice , the special energy 
value E = €a is determined by considering the AA cluster alone. Many 
other Id lattices containing these AA dimers may have an extended state 
at this same energy. To illustrate this let us consider a periodic lattice of A 
atoms. Here the AA correlation is trivially present, and therefore this system 
has an extended state at the same energy. However, while all eigenstates at 
energies other than ca are localised in the random dimer case, here we have 
a whole band of extended states, the existence of which cannot be inferred 
by considering the AA dimers. Although Bloch's theorem enables us to find 
this band directly, we may alternatively think of determining this band by 
considering higher order clusters consisting of triplets, quadruplets . . . etc. 
of A atoms, as has been discussed by us earlier 5 . There is no advantage in 
adopting this point of view in this elementary example. However, in aperiodic 
chains Bloch's theorem does not apply, and one has to use this idea in order 
to determine the bands or mini-bands of extended states. 

The clusters contain only partial information about the entire lattice, and 
consequently, we have the same set of energy levels for systems all of which 
contain a particular cluster but which have different long-range compositions. 
The detailed differences in the electronic spectra in these systems are due to 
correlations at length scales beyond the cluster size. In order to bring out 
these differences we have to examine correlations among larger and larger 
blocks of atoms in the chain, by looking at the systems at increasing length 
scales. In this process the long range features of the lattice get gradually 
included and we end up by obtaining the entire spectrum of extended states. 
The real space renormalization group (RSRG) method is particularly suited 



for studying this kind of problem. 

In this paper we consider a period-doubling (PD) chain 10 ' 12 as a prototype 
example to illustrate the above clustering ideas. These ideas are quite general 
and are also applicable to other aperiodic systems. A PD chain, to our 
mind, is a very interesting candidate for investigation beacuse the trace map 
associated with this chain leads to a polynomial invariant which is likely to 
ensure the critical nature of all the eigenstates as well as a Cantor set energy 
spectrum. As we will show in this paper, inspite of the existence of such an 
invariant, this lattice can sustain an infinity of extended eigenstates which 
coexist with the otherwise Cantor set nature of the energy spectrum. 

The few analytical studies that are available on the PD lattice have all 
been made with a simple on-site tight-binding version of the lattice hamil- 
tonian. Bellissard et al. 11 have shown that the spectrum of a tight binding 
Schrodinger equation in which the potentials are distributed following a PD 
sequence is purely singular continuous and supported on a Cantor set of zero 
Lebesgue measure. In an earlier article Severin et al. 12 analytically studied 
the case of a special class of quasiperiodic sequences generated by the sub- 
stitution rule Sl+i = Sl-iSl-iSl which yielded the standard copper mean 
lattice and variations of PD lattices on tuning the initial choice of the basic 
building blocks. They reported the existence of periodic eigenstates on such 
lattices at some special values of energy. From our point of view, the work of 
Severin et al. is particularly interesting, as they give an analytical method 
for locating the eigenvalues and the corresponding extended eigenstates for 
the PD sequence. The validity of their approach, however, is strictly limited 
to the use of the on-site model. Apart from this, their analysis crucially 
depends on a special choice of the value of the pseudo-invariant 12 7 for PD 
chain equal to 2, and it is not clear how this analysis can be extended to 
other values of 7. 

Interestingly, in all the works on the PD chain reported so far, the clus- 
tering effect that we have discussed above, although present, has neither 
received adequate emphasis, nor has been made the focal point for analysing 
the nature of the eigenf unctions. As will be seen, the clustering idea enables 
us to reduce the problem of identification of the nature of eigenfunctions 
on any quasi-periodic lattice to its essential simplicity. Thus, for instance, 
by merely looking at the Fibonacci chain it is possible to infer that there 
are no extended states in general, without going into any further analytical 
consideration, since there is no clustering in this chain 8 . 



One of the aims of this paper is to investigate the nature and character 
of the electronic eigenfunctions with more general hamiltonians. Such an ex- 
tension is necessary because any experimental realisation of a quasi-periodic 
chain will not be restricted to on-site variations of parameters only. As a first 
step we therefore look at a tight-binding model hamiltonian for the PD chain 
which includes simultaneous variations in both on-site and hopping matrix 
elements. This increase in the level of complexity is not readily amenable 
to analysis through the transfer matrix approach alone. For the PD case we 
find interesting cross-over behaviour in the nature of the eigenstates from an 
all states critical picture to a situation where we have an infinity of extended 
states, depending on the region of parameter space in which we work. We 
next consider the problem of the motion of an electron in an array of 5- 
function potentials whose strengths are distributed in accordance with the 
underlying quasi-periodicity of the lattice. This is a multi-band problem, 
and we find that the clustering idea still applies. We may regard this case as 
approximating the experimental situation more closely than the single band 
model based on the tight binding approach. 

We organise this paper as follows. In the next section we give a somewhat 
detailed exposition of the clustering idea in the context of general single band 
tight binding models of a PD chain. We also present numerical results for 
the extended eigenfunctions obtained on the basis of the present theory. The 
extension to the multiband case is given in section III. Section IV deals with 
the multifractal analysis of the wave functions. 

II. The Tight Binding Model 

A portion of the period doubling chain is shown in Fig. 1(a). The sequence 
in which the long (L) and short (S) bonds are arranged in this chain is 
obtained by successively using the substitution rule 10 L — > LS and S —* LL. 
For describing the electron states in this lattice we use the single band tight 
binding hamiltonian 

H = J2ei\i><i\+12 t ij\ i>< J\- W 

Since we are interested in locating clustering effects at different length 
scales, it is essential to label the sites and the hopping matrix elements ap- 
propriately. This is necessary, because in order to implement the RSRG 
decimation scheme, the self-similarity inherent in the lattice has to be pre- 



served at every stage of renormalization. Such labelling was discussed by 
us earlier for the Fibonacci lattice and the copper mean lattice 5 ' 8 . There is 
a basic difference in the renormalization procedure followed in the transfer 
matrix method and that in the present approach. Our approach is based 
on the global decimation procedure discussed by Southern et al. 14 , in which 
we consider an infinitely long chain, and decimate a subset of sites so as to 
produce a self-similar chain on an inflated length scale. We find that we 
have to assign at least three different labels to the site energies and two to 
the hopping matrix elements to make this decimation possible in the PD 
chain. Consequently we identify three types of sites, which we label by a, /3 
and 7 corresponding respectively to sites between two consecutive long (L) 
bonds, a long (L) bond and a short (S) bond, and a short (S) bond and a 
long (L) bond. The site energies in Eq.(l) therefore assume the values e a , 
e@ and e 7 respectively. The hopping matrix elements across the long bond 
(L) is taken as t^ and that across the short bond (S) as tg- As we shall see, 
the site energies and the hopping matrix elements of the renormalized chain 
at any stage carry these very five labels and no more. Of course later we 
shall discuss the very interesting case of a still more general model which has 
four types of hopping matrix elements across the long (L) bond depending 
on the vertices connected, and whose behaviour is quite different requiring a 
separate treatment. 

Let us begin with the case where we have three site energies e a , ep and 
e 7 and two hopping matrix elements ti and ts- The eigenfunctions for such 
a lattice may be calculated by the standard transfer-matrix method. The 
Schrodinger equation for the hamiltonian (1) can be cast in the form <p n+ \ = 
M n cf) n , where 




*n,n— 1 



and M, 



1 1 




Here ip n denotes the amplitude of the wave function at the n-th site and M n 
is a 2 x 2 transfer matrix. In the PD chain we have three different kinds of 
transfer matrices M a , Mg and M 7 , and n is related to 0o by a product of 
these three matrices following the PD sequence. 

By inspecting the PD chain (Fig. 1(a)) we see that the a-sites always occur 
in pairs whereas, the (3 and 7 sites always occur as a doublet f3j. Thus for 
this lattice the pair of sites a — a constitutes a cluster in the sense discussed in 



the Introduction. For this chain the string of transfer matrices typically looks 
like. . . M M .Ml.M M .M M .M M .Ml . . ., where M M = M y Mp. We notice that 
the matrix M a is unimodular. The pure on-site model results as a special 
case by taking e a = e 7 = €a, tp = £b an d ti, = t s while the transfer model is 
obtained by taking e a = ep = e 7 = e and t L ^ t s . 

Let us first discuss the general case of a one dimensional chain which 
contains m-component clusters, that is clusters each containing m identical 
building blocks. Let the transfer matrix for every building block be a 2 x 2 
unimodular matrix denoted by M. The composite transfer matrix for the 
cluster is then M m . By a straightforward application of the Cayley-Hamilton 
theorem for the ra-th power of a 2 x 2 unimodular matrix M it can be shown 
that 

M m = C/ m _i(x)M - U m . 2 (x)I (2) 

where, x = (l/2)TrM. U m (x) = sin(m + 1)6?/ sin #, with 9 = cos^rr is the 
m-th Chebyshev polynomial of the second kind 5 . For a value of the energy E 
for which U m -i(x) becomes zero, M m reduces to — f/ m _ 2 (x)J, that is to say, 
the transfer matrix for this cluster behaves essentially as the identity matrix 
at this energy. Thus at this energy the entire lattice does not feel the presence 
of the cluster defined by the transfer matrix M m . If the remainder of the 
lattice forms a periodic chain, there will be extended states at this energy 
provided, this energy is an allowed one. For allowed states wave functions 
do not diverge at infinity. 

For the PD chain M = M a and m = 2 and therefore setting U\(x) = 
we obtain E — e a . For this energy the whole lattice effectively behaves as 
a periodic diatomic linear chain with each unit cell containing a (3 and a 7 
atom. If the energy E = e a happens to be within the allowed band of this 
diatomic lattice, then this energy is an allowed one. The condition for this 
is (l/2)|TrM^ 7 | < 1. 

We now proceed to determine other energy eigenvalues which may give 
rise to extended eigenfunctions for the entire chain. From the Schrodinger 
equation for the PD chain (see Fig. 1(a)) we obtain the following hierarchy of 
equations for the amplitudes of the wavefunction 



(E - e 7 )V'-i = t L % + tsi>-2 
(E - e a )ij Q = tLi>i + t L ip-i 



(E-eJip! = t L ^ 2 + t L iJo (3) 

(E - €p)i) 2 = t 5 ^ 3 + t L lf)i 

(E - e 7 )^ 3 = t L i) 4 + t s ip 2 



If we generate a renormalized PD chain by decimating an appropriate set 
of sites ( see Fig. 1(a) ),then we obtain a self-similar hierarchy of equations, 
but with renormalized parameters corresponding to the inflated chain shown 
in Fig. 1(a). The topology of the chain is preserved as a result of this trans- 
formation and the corresponding scale factor is two. The renormalized site 
energies and the hopping integrals are found to be 

/ 9 9 

ep = e 7 + LO a t L + Wpt s 

e 7 = e a + (u a + up)t 2 L (4) 

ti — UptLts 

where u>i = l/(E — ei), i = a, (3, 7. Corresponding to this transformation it 
can be shown that there exists the following polynomial invariant 

J " W~s + ' (5) 

whose value remains unchanged under the transformation Eq.(4), as can 
be easily verified. Since the renormalized chain is still a PD sequence, we 
again find a-sites occuring as a — a pairs and /3 — 7 sites forming (3 — 7 
doublets. This means that the a — a clustering effects are also present on 
this inflated length scale. With respect to this renormalized lattice there are 
extended eigenstates at E — e' a (when M 2 = —I), provided for this energy 
k' = (l/2)|TrM / g 7 | < 1, a, /3 and 7 now referring to the renormalized chain. 
Since the a — a clustering is present at every length scale, one will find, upon 
repeated renormalization, a number of extended state-energy eigenvalues by 
solving the equation E = e£> , and by checking that the roots of the equation 
satisfy the condition k^ < 1, where the trace k^ is evaluated with the 
renormalized parameters at the n-th stage. Interestingly, for the PD chain 
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we find that the value of k^ becomes equal to X — 1 for energies which are 
the roots of the equation E — e^ = 0, where, X is the invariant defined in 
Eq.(5). This implies that k^ is also independent of the generation index n 
for these special energy values. Therefore, if the initial parameters are such 
that k^ is greater than one, then there will be no extended states at any level 
of renormalization, and the roots of the equation E — e^ = will always be 
in the gaps of the entire band. On the other hand, if for the initial choice of 
parameters, /t*- -* is less than or equal to unity, then all the subsequent «^'s 
will also be < 1, and thus all the solutions of the equation E — e^(E) = 
will be allowed ones. 

It is to be appreciated that the a — a clustering at larger length scales 
amounts to including, so to say, the effects of larger and larger segments of the 
original chain into the a-subclusters. The a— a correlation at all length scales 
can only be revealed by renormalization group methods. While in the original 
lattice this correlation is directly visible (Fig. 1(a)), higher order correlations 
due to a — a pairing imply underlying complex correlations between atoms, 
which is not apparent from mere inspection of the original lattice. 

With this background we are now in a position to discuss several mod- 
els of the PD chain which are obtained by assigning different values to the 
hamiltonian parameters. Specifically, we shall discuss the following models: 

(i) The on-site model 

In this model the hopping matrix elements are all taken to be equal and 
the site energies are of two types, e^ and cb, arranged on a lattice following 
the PD sequence. Severin et al. 12 in their analysis showed that this chain 
supports extended states which display periodicity of periods 4, 8 . . . etc. in 
units of lattice spacing. They made use of the fact that a pseudo-invariant 
associated with this model has a value equal to 2, which is then used to 
establish the periodic nature of the solutions. From our point of view, we 
recover the on-site model by putting e a — e 7 = ca, eg = e_e and th — ts — t. 
For numerical calculations we choose 6 a = 1, 6b — — 1, t — 1 in suitable 
units, which are the same as those used by Severin et al., in order to afford 
a comparison with their calculation. At the basic level the extended state 
occurs at E — e^ — 1, and the wavefunction for this energy is found by the 
transfer matrix procedure beginning with the value ipi = 1 and ip = for 
a 7th generation PD chain with 128 atoms. The wavefunction is shown in 
Fig. 2 (a) and is identical with that of Severin et al. 12 having a periodicity of 



4 units of lattice spacing. This energy value yields for the quantity k^ the 
value -1, corresponding to the edge of the unperturbed /?7 band of energies. 
Since this quantity equals X— 1, the trace values k^\ kP\. . . k^ . . . , are all 
equal to -1. If we now consider the renormalized chain, then the extended 
states occur at energies obtained from the equation E — e' a (E) = 0, the roots 
being ±a/3. The wavefunction for E = \^3 is shown in Fig.2(b), with a 
period of 8 units, and is again identical with that of Severin et al. 12 At the 
next level we have the energies by solving E — e'^{E) = 0, which yields the 
roots E = -1.82595501, -0.15244466, 1.38553731 and 2.59286237. In Fig.2(c) 
we give the wavefunction for E = 2.59286237. This function is seen to have 
a period of 16 units. Finally, in Fig.2(d) we plot the wavefunction for the 
energy E = 2.60380559, obtained from the next level of renormalization. 

As was pointed out by Severin et al., the underlying reason for the ob- 
served periodicity of say, the eigenfunction corresponding to E — 1 having a 
period of four lattice spacings, is related to the vanishing of the amplitude at 
every fourth site on the chain. Alternatively, we note that at this energy, the 
hierarchy of Eqs.(3) are indistinguishable from that of an infinite periodic 
lattice of alternate A and B atoms. Similarly, the equations determining the 
eigenfunctions of period 8 are identical with that of an ordered lattice with 
unit cell consisting of ABAA atoms. In the same manner it is found that 
every periodic solution in the on-site model is identical with the solution of 
some periodic chain with a suitable unit cell. Thus it turns out that in this 
case we may associate a whole array of periodic lattices with different unit 
cells with the extended eigenfunctions, each periodic lattice being in one to 
one correspondence with one of the eigenvalues of the PD chain. A similar 
observation has been made recently by Oh and Lee 16 . 

(ii) The transfer model 

As mentioned before, the transfer model is obtained by choosing the start- 
ing values e a = ep = e 7 , and t^ ^ t s . For numerical work we choose the 
values e a = e/3 = e 7 = 0, t^ = 1 and t s = 2. For this case the energy for 
the extended state E = lies in the central gap of the band, with a value 
of k^ = —1.25. Again, from the invariance of the trace under the RSRG 
transformation, it follows that k^ = —1.25 for all n, that is, all the energy 
values are disallowed and thus there are no extended states in this model. 
In fact, in this model the value of k^ for any choice of ti and ts is of the 
form «( n ) = (1/2) (x+ l/x), where, x = t L /ts, and thus k^ is always greater 
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than unity. Since the value of X in this case is non-zero, the eigenstates are 
all critical supported on a Cantor set of zero measure. 

(iii) The mixed model 

Let us begin with the case where all the site energies are unequal, and 
there are two different hopping matrix elements £l and ts- The transfer ma- 
trix M a for this case is unimodular, and hence the energies for the extended 
states are again obtained from the equation E — e^(E) = 0. The band of 
energies allowed by the periodic f3^f chain defined by the initial parameters 
lies in the range 

'ep + e 7 + y/fo - e 7 ) 2 + A(t s + t L ) 2 e p + e 7 + ^(e - e 7 ) 2 + A(t s - t L )< 




(ep - e 7 ) 2 + 4(t s - t L f e p + e 7 - J(e p - e 7 ) 2 + 4(t s + t L ) 



\2 



If the value of e a lies outside these energy intervals, then k^ ' becomes greater 
than unity at the energy E = e a and there is no extended state at any level of 
renormalization. With the progress of renormalization the periodic /?7 chains 
with renormalized parameters will yield more and more fragmented bands, 
and the energy values obtained as solutions of E — e^ (E) = will always be 
in the gaps of these fragmented spectrum. Since the value of X is non-zero, 
all the states will be critical even in the presence of a — a clustering. On 
the other hand, if the value of e a lies inside one of the initial energy intervals 
defined above, we always have k^ < 1, and the states will be extended. The 
extended character of each wavefunction has been checked by a mutifractal 
analysis (see section IV). There is thus a cross-over in the behaviour of the 
eigenstates depending on the choice of the initial parameters. Confining 
our attention to the regime of extended behaviour, we find an interesting 
systematics in the nature of the eigenfunctions as we consider the roots of 
E — e^\E) = for increasing values of n. Choosing e a = 4, ep = and 
e 7 = 3, tL = 1 and ts = 2, we find that all the wavefunctions are non-periodic, 
as can be seen from Fig. (3), where we have displayed the eigenfunctions at 
the energies 4, -1.19258240, 5.45570651 and 3.56527726 arising from the 
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first four levels of renormalization. In Fig. (4) eigenfunctions are shown for 
energies -0.60581006, 4.14562305, 4.65276530 and 5.45820519, which arise 
from the solution of E — e^ (E) = 0. As may be seen, the amplitudes of these 
non-periodic functions do not decay as we go from one end of the chain to the 
other. We have tested the non-decaying character of the wavefunctions for 
chain lengths upto 2 18 atoms, although we have presented the amplitudes for 
much smaller chain lengths for convenience. A very interesting feature that 
may be noted is that eigenfunctions corresponding to neighbouring energies 
possess, in general, entirely different profiles. This feature is contrary to 
what one obtains in a periodic lattice, and one may regard this behaviour as 
a manifestation of the quasiperiodic character of the lattice. 

We now go on to the discussion of a still more general model, namely 
one in which, in addition to having three different site-energies, we ascribe 
four different values to the hopping integral across the long bonds connecting 
aa, a[3, 7a and 7/3 pairs of sites. While the primary reason for considering 
this model is that it represents the realistic situation most closely within the 
framework of the single band nearest neighbour tight-binding hamiltonian, an 
additional reason is that even though a — a clusters are seen to be present, 
the matrix M a will have two different values depending on the immediate 
neighbours of an a site. Neither of these two M a matrices is unimodular, 
and therefore one cannot apply the result of Eq.(2) to this case. Inspite of 
this, the present model gives rise to an infinity of extended states. To see this, 
if we consider only the subset of a sites in this PD chain, then these sites are 
found to form a lattice in which two different bond lengths are arranged in 
the copper- mean sequence 5 . In the same manner, each of the subset of /3 and 
7 sites forms a copper-mean sequence. If we now decimate the original PD 
chain so as to retain only the j3 sites (the same could have been done for the 
a or 7 sites), then we may relabel the sites of the resulting lattice following 
the prescription of the copper- mean chain discussed by us earlier 5 , so that 
the original PD chain transforms into a copper mean lattice (see Fig. 1(b)) 
with four site-energies e^ M , e? M , e^ M and ef M and two different hopping 
integrals tjj M and t$ defined as 

e« M = e/j + st|w 7 + qt 2 a pu a 

CM _ 1 .2 1 .2 

e /3 — e /3 + ts^i + qt Q/3 uj a 

e c ' M = ep + st^ + t'u, (6) 
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C CM _ j_ .2 , , , - 2 , , 

tg = ty a tgUJ^i 

where oVs have been defined previously, p, q, r and s are given by 

V = (l-t^o^r 1 

Q = i 1 ~ Ptla^l)' 1 

r = (1-iL^r 1 
s = {1 - rt^UaU^Y 1 . 

Very interestingly we find that as a result of this transformation the 
hopping matrix elements for the effective copper-mean chain takes only two 
different values depending on the two distinct bonds in this renormalized 
lattice, even though we had started with the most general version of the 
PD chain. In this renormalized copper-mean lattice we have again the a — 
a clustering as discussed in Ref.5, and this leads to a whole hierarchy of 
extended states, provided (l/2)|Tr(M 7 M (5 M /3 )| < 1. 

Two comments are in order at this stage. Firstly, it is impossible to 
discern this clustering effect at the level of the original PD chain, and the 
full spectrum of extended states can be obtained only through this initial 
transformation to a copper-mean lattice. This means, that the original PD 
chain does implicitly contain a complex clustering involving all the a, (3 and 
7 sites, but which does not show up at the initial stage. Secondly, it may 
seem that one could have straightaway applied this transformation to any of 
the models discussed earlier. However, it should be appreciated that there 
is a basic difference between the RSRG transformation for a PD lattice and 
that for a copper-mean chain. While there exists a polynomial invariant X 
associated with the PD chain, there is no such invariant for the copper-mean 
case. The value of X for the most general model is found to be 

(E(t aa - t 7/3 ) - t aa e 7 + e a t jfS )(e a - ep) - t s t aa (t s - 2t 7/3 ) - t aa t lf3 

(7) 

In determining X we have set the initial hopping integrals t aa — t a p and 

t ia = t 1 @. This is, of course, not a restrictive condition, since starting with 

the most general model with four different values for the hopping integral 
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for the long bonds, we get, upon decimating alternate sites following the PD 
renormalization rule, only two different values for this hopping ,which we 
may choose to be t aa and t 1 p. The quantity X is a linear function of energy, 
in contrast to the any of the other models discussed earlier. By setting 1=0 
(see Ref.7 and 8), we immediately find that there is an additional extended 
state at an energy 
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We have computed this wavefunction numerically with e a — 1, ep — 0, 
e 7 = 0.1, t aa = 1, t ya = 1.1 and ts = 2, for which E is found to have a value 
1.428571428. This extended state cannot be detected from a consideration 
of the effective copper mean lattice. The character of this extended state 
is entirely different from the other extended states in the PD lattice, if one 
looks at the pattern of flow of the hamiltonian parameters in the copper-mean 
chain with renormalization. As has been pointed out by us earlier 5 , for an 
energy obtained from the equation E — e^ (E) = for the effective copper- 
mean lattice, the values of ep and e 7 become equal after a certain number of 
iterations depending on the value of n, and remain so in subsequent iterations. 
Moreover, the hopping integrals do not flow to zero, which is an evidence of 
the extended nature of these eigenfunctions. On the other hand, for the 
special energy value defined in Eq.(8), ep is never found to become equal 
to e 7 at any stage of iteration, although the hopping integrals do remain 
finite and non-zero, indicating that this state is also extended. The chaotic 
behaviour of the flow pattern may be said to characterise this eigenstate. The 
corresponding wavefunction is shown in Fig. (5), and its extended character 
is confirmed through the multifractal analysis given in section (IV). 

III. The Multi-Band Model 

Let us now consider the continuum version of the Scrodinger equation in 
one dimension with a quasi-periodic potential V(x). We may transform this 
equation with any arbitrary potential to a discrete set of difference equations 
using the standard Poincare map 13 . If we label the points on a Id lattice by 
the integers Xi, then we can find a recursion relation for the amplitudes of 
the wavefunction at three consecutive points of the above set. In the usual 
language of tight binding theory, the "hopping matrix element" connecting 
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amplitudes on neighbouring points Xi can be expressed in terms of the Wron- 
skian matrix corresponding to the solutions of the Schrodinger equation in 
that interval. We omit the expressions here and refer to the literature for 
details 13 . For simplicity, let us apply this theory to the case of an array of 
(^-function potentials sitting at the sites % with strengths \ distributed ac- 
cording to the PD sequence. The Poincare map leads us to a set of difference 
equations of the following form connecting the amplitudes at three successive 
points: 

ipn+i + VVi-i = [2 cos q + (Xn/q) sin q}^) n (9) 

where q = J(E). Here \ n gives the strength of the potential at the n-th site. 
The energy eigenvalues forming bands of extended states correspond to the 
values of the wave vector q satisfying the inequality 

2cosg+ (A n /g)sing < 2 . (10) 

In order to extract information about the spectral nature as well as the char- 
acter of the eigenfunctions in this case we can recast Eq.(9) in the following 
form : 

(E' - e n )i) n = V>n-1 + i>n+l , (11) 

with E' = 2cosg and e n = —(X n sinq)/q. This is now in the familiar form 
of the single band tight binding hamiltonian where E' now plays the role 
of the 'energy' and e n can be interpreted as the 'effective' on site potential 
distributed according to the PD sequence. We may now assign three different 
labels a, (3 and 7 to the sites of the chain and adopt the same procedure as 
described earlier for finding the extended eigenstates. Setting E' = e^\ 
where, n implies the n-th stage of renormalization, we have been able to 
locate the eigenvalues corresponding to the extended eigenfunctions for the 
multiband PD chain. There is however, an important aspect to note, viz. 
setting E' = e a even at the initial stage we arrive at an equation of the form 

2/\ a = -(l/q)tanq. (12) 

Now tang is a n -periodic function and it takes all values in the interval —00 
to 00 for any value of A Q . Therefore, at each stage of renormalization the 
solution of the equation E' = t^ will give rise to an infinite countable set of 
energy eigenvalues for which the PD chain will be totally transmitting. That 
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is to say, we have infinite number of extended eigenf unctions at all scales of 
length. This result is to be contrasted with the case of a tight binding random 
dimer model which exhibits a resonance state at a unique energy, or with the 
results for the copper-mean chain or the PD chains in the tight binding 
versions where, at each stage of renormalization one generally enumerates 
a finite number of extended eigenstates. For numerical calculation of the 
wavefunctions we have taken \ a = A 7 = 1 and Xp — — 1. At the very 
initial stage, by setting E' = e a we obtain an infinite number of allowed q 
values. We have explicitly calculated the wavefunction for q = 1.836597203 
and the result is shown in Fig.6(a). This wavefunction is real and periodic 
with a period four. The amplitudes follow the sequence 1,0,-1,0 for the 
first four sites and this pattern is repeated for the rest of the lattice. The 
discontinuity of the derivative of the wavefunction at lattice points x = 3na, 
n = 1,2,3... have been shown encircled in Fig.6(a). From the next level of 
renormalization, we have chosen q = 2.37738316 from the set of allowed q- 
values, and have plotted the wavefunction, which has a periodicity of eight, 
in Fig.6(b). The wavefunctions for other levels of renormalization may be 
obtained in the same manner. 

IV. Multifractal Analysis of the Wavefunctions 

All the wavefunctions given in this paper have been subjected to a mul- 
tifractal analysis to verify their extended nature. Following the standard 
procedure, a "partition function" is defined as 

Z(Q) = J:\H Q (13) 

1=1 

where N is the number of sample points. For sufficiently large N, Z(Q) may 
be expressed as 

Z(Q) = e r(Q) (14) 

where e is the interval 1/N, and t(Q) is an index from which the multifractal 
index a and the corresponding fractal dimension /(a) are obtained as 

a = ^?1 and f{a) = Qa - r{Q) . (15) 

As is known 15 , the standard algorithm is however not of much practical 
utility unless the sample size is inordinately large. With all finite, reasonably 
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large, sample sizes one always obtains a distribution in the values of a and 
f(a). This is a finite-size effect and the converegence to the true behaviour 
occurs only logarithmically slowly with increase in the sample size 15 . Thus 
the standard procedure cannot be regarded as a reliable method to analyse 
finite size samples. As has been explicitly shown by Godreche and Luck 15 , 
a much better method to test the extended character is to look at the be- 
haviour of the curvature of a — f{a) curve at Q — 0. As has been shown by 
them, the quantity |1//"(o0|q=o, which is a measure of the curvature at the 
maximum of the a — f{ot) graph, should diverege linearly with log(iV) for 
extended states, N being the sample size. On the contrary, if the given set 
is a true multifractal, the curvature should gradually saturate to single value 
with log(iV). In our case we have used this prescription to test our wavefunc- 
tions, with a maximum sample size of 2 17 values of the amplitudes. In every 
case, we observe the predicted linearity in the dependence of |1// // (ok)|q=o 
with log(iV). In Fig. 7 we give the results for three of the wavefunctions, 
including the 'chaotic' wavefunction shown in Fig.5, the graphs for all the 
other wavefunctions being entirely similar. Our analysis, therefore provides 
confirmatory evidence that the states obtained on the basis of the present 
work are true extended states. 
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FIGURE CAPTION 

Fig.l (a) A section of the period-doubling chain. The a, (3 and 7 sites are 
indicated respectively by a full circle, a square and an open circle. 
a — a clusters are shown encircled by dotted lines. The renormalized 
chain under a period-doubling decimation is also shown, (b) Effective 
copper-mean chain obtained by eliminating the a and 7 sites in the 
period-doubling chain. The a, /3, 7 and S sites are marked by a full 
circle, square, open circle and a triangle in the copper-mean chain. 
a — a clusters in the copper- mean chain are also indicated. 

Fig. 2 Wavefunctions for the period-doubling chain in the on-site model for 
e a = e 7 = 1, ep = —1, ti, = 1 and t s = 2. For (a) E — 1, (b) E — \/3, 
(c) E = 2.59286237 and (d) E = 2.60380559. All energies are in units 
oit L . 

Fig. 3 Wavefunctions for the period-doubling chain in the mixed model for 
e a = 4, ep = 0, e 7 = 3, t L = 1 and ts = 2. For (a) E = 4, (b) 
E = -1.19258240, (c) E = 5.45570651 and (d) E = 3.56527726. All 
energies are in units of t^. 

Fig. 4 Wavefunctions for the period-doubling chain in the mixed model. Pa- 
rameters are the same as for Fig. (3). The energies are for (a) E = 
-0.60581006, (b) E = 4.14562305, (c) E = 4.65276530 and (d) E = 
5.45820519. 

Fig. 5 Wavefunction for the general model : e a = 1, ep = 0, e 7 = 0.1, t aa = 
t a p = 1 and t ia — t^p = 1-1, and ts = 2, for E = 1.428571428 obtained 
from Eq.(8). The energies are in units of t aa . 

Fig. 6 Wavefunctions for the period-doubling chain in the multi-band case for 
\ a = A 7 = 1, A^ = —1 for (a) q = 1.836597203 corresponding to a 
period four and (b) q = 2.377380316 for a period eight. The energies 
are measured in units such that h /2m = 1. The discontinuities in 
ip'{x) are shown encircled. 

Fig. 7 Plot of |l//"( a )lQ=o (curvature) against log(iV) corresponding to wave- 
functions in Fig. 2(d) , Fig. 3(d) and Fig. 5 are shown in (a), (b) and (c) 
respectively. The maximum sample size has been taken as 2 17 . 
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